Some applications of the Faddeev-Yakubovsky equations to the cold-atom physics 
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We present some recent applications of the Faddeev-Yakubovsky equations in describing atomic 
bound and scattering problems. We consider the scattering of a charged particle X by atomic 
hydrogen with special interest in X = p, e^ , systems of cold bosonic molecules and the bound and 
CO ' scattering properties of N—3 and N=4 atomic ''ife multimers. 
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^ ■ I. INTRODUCTION 

^ ; The Faddeev-Yakubovsky (FY) equations constitute a rigorous mathematical formulation of the quantum mechan- 
ical N-body problem in the framework of non relativistic dynamics. They allow obtaining exact solutions - in the 
numerical sense - of the Schrodinger equation for bound and scattering states, in principle for an arbitrary number of 
particles. They were first formulated as integral equations in momentum space by Faddeev in the early sixties [1], in 
I the context of 3-nucleon problem with short range interactions, and were generalized some years latter to an arbitrary 

I number (N) of particles by Yakubovsky . 
^ . A substantial progress in their numerical solution was made when the boundary conditions in the configuration 
O ■ space were established for N=3 in and for N>3 in Q. A reformulation of the original Faddeev equations, 

I ^ I , allowing to incorporate long range Coulomb like interactions was derived in 0, Q • A brief review of the different steps 
in elaborating the Faddeev-Yakubovsky approach to N-body problem can be found in . 
I Until now only N=3 and N=4 systems have been explored, although there is no any formal nor practical impediment 

^ to go beyond. Nowadays very accurate FY solutions have been obtained for the bound state problems on hadronic, 
nuclear and atomic systems, in configuration as well as in momentum space [1343] i still with a smaller accuracy than 
I the variational methods. The real advantage in using the FY equations, which was also their original motivation, 
is in describing on the same foot the scattering states. They play there an unavoidable role, especially when many 
channels are open including the break-up ones (dissociation into more than two clusters). At present, there exist 
J;^^ • three-body reliable scattering results in nuclear and atomic physics including the break-up process [T5l4T9| . In the 
four-body case, the progress remains limited to elastic and rearrangement 1-1-3 and 2-1-2 particle channel in nuclear 
and atomic problems [ill, Us, 20-23]. 

We present in this contribution some recent applications of the Faddeev-Yakubovsky formalism to describe cold 
atomic structures. A short introduction to the Faddeev equations is given in Section |lll Section UTT] is devoted to 
present some results of a charged particle interacting with atomic hydrogen, with special interest in the p — H (or 
, molecular ion) and — H systems. We study in Section ITVl universal properties of the weakly bound N=3 and 

C^_- N=4 cold bosonic molecules. Some results concerning "'He atomic multimers are given in Section |Vl The last section 
is devoted to draw some conclusion and perspectives. 
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II. THE FORMALISM 



We will develop in what follows the main ideas of the Faddeev equations for solving the three-body problem in 
configuration space. This form is the best adapted to deal with atomic problems. Let us consider three different 
spinless particles with masses rrii, coordinates and denote by {xi,yi) the three different sets of Jacobi coordinates 
defined by 



Vi = 



mo{mj + ruk) 

l 2nii{mj + rrik) 
irioM 



i = 1,2,3 



rrijrj + mkTk 



rUk 



(1) 



(2) 



where {ijk) denotes a cyclic permutation of the particle numbers (123), M — mi + m2 + is the total mass of the 
system and mg an arbitrary mass to be fixed at convenience. We will also assume that these particles interact via 
pairwise short range local potentials Vi{xi) vanishing at some distance Ri. 





FIG. 1: Sets of three-particle Jacobi coordinates. 

The dynamics of the system can be alternatively described by {^^1^=1,2,3 or by choosing one of the Jacobi coordinate 
sets (xi, i/i), completed by the center of mass coordinate R 

MR — rriifi + 7712^2 + "1,3^3 
The three-body Hamiltonian is assumed to have the form 

n = no + v (3) 

where Ho is the kinetic energy which can be written in one of these forms 



Ho 



and V the total potential 
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In order to remove the center of mass motion, the choice of a Jacobi set is mandatory. Only in terms of it, the 
three-body wavefunction 5', eigenstate of the total three-body Hamiltonian ([3]), factorizes into 

where ^{xi,yi) is an eigenstate of the intrinsic Hamiltonian 



H 



mo 



(5) 



the only one that we are going to consider from now. 

As it can be already seen from none of the Jacobi sets is neither privileged nor fully satisfactory. The "non 
interacting" region of particles 2 and 3 is simply given by xi > Ri but it is difficult to define this region in terms of 
say (x2,y2). 

On the other hand, if we wish to describe for instance the scattering of particle 1 on a bound state of particles 2 
and 3 - denoted symbolically by 1(2,3) - it will be natural to chose the coordinate set {xi,yi). However the final 
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state can contain, together with the initial state, a superposition of channels 2(1,3) and/or 3(1,2) which can be hardly 
described in terms of the same coordinate set. For the bound states and break-up channels the three ensembles will 
appear naturally on the same footing. 

The Faddeev equations are based on a partition of the total wave function on as many components as two-body 
asymptotic channels: 

<I> = $1 -)- <i>2 -f <f>3 

It is straightforward to see that the three-body Schrodinger equation ([5]) is equivalent to the set of coupled partial 
differential equations for the Faddeev components 

[E - Ho-Viixi)]^i{xi,yi) = V"i(a;i)[*2(^2,y2) + $3(^3,^3)] 
[E - Ho -V2ix2)]'^2{x2,y2) = V'2(a;2) [$3(^3,2/3) + $1(^1,2/1)] 

[E-Ho-V3{x3)]^3{x3,y3) - 1/3(2:3) [*i(2?i,yi) + $2(2^2,^2)] (6) 

The coupling is ensured by the right hand side. It is highly non local due to the linear relations between two different 
sets of Jacobi coordinates: Xa^xp, yp), ijaixp, yp). In the non interacting region, Vi — 0, the three Faddeev equations 
decouple and the boundary conditions for each component take a simple form when expressed in their own Jacobi 
coordinate set. 

These boundary conditions are more easily implemented in terms of the reduced Faddeev components 4>i defined 
by: 

(f>i = Xiy.^i (7) 

and take the following Dirichlet form: 

- they vanish for Xi — and yi — 0: 

(j>i{xi = 0,yi) = (j)i{xi,yi = 0) = (8) 

- for a 3-body bound state they decrease exponentially in all the directions. In practice one can force them to 
vanish at sufficiently large distances 

4^i{Xi 2^ "^niaxj yi ^ ymax) — (9) 

- for an open i-l-(jk) elastic or rearrangement scattering, the i-th Faddeev component splits into the product of 
the two-body bound state wave function ipi{xi) of the particle pair (jk) and the scattering wave of particle i 
with respect to the center of mass of this pair XiiVi)'- 

(l)i{xt,yi) = (pi{xi)xi{yl) (10) 



for the break-up reactions at large values of the hyperradius p — \/ xf + yf one overimposes to (llOp the behavior 



Mxi,yi)^-^ (11) 



As mentioned in the Introduction, the original Faddeev equations, above formulated, are not suitable for the 
Coulomb scattering problems. The reason is that the right hand sides of eq. ^ do not decrease fast enough to ensure 
the decoupling of Faddeev amplitudes in the asymptotic region and to allow unambiguous implementation of the 
boundary conditions. In order to circumvent this problem, Merkuriev 0, |1| proposed to split the Coulomb potential 
V into two parts, V = V + V'' , hy means of some arbitrary cut-off function rj: 

V'{x,y)^V{xUx,y) (x^y) ^ V{x)[l - r^{x,y)] (12) 

One is then left with a system of equivalent equations 

[E-Ho~W^- y/] = + $fe] = + V,' + V/: (13) 

with the right hand side containing only short range contributions (K,) and with a 3-body potential (Wq,) in the left 
hand side to account for the long range parts. As it will be demonstrated in the following, this approach was found 
to be very efficient in calculating the positron-positronium (Ps) and positron- hydrogen (H) cross sections [TtI . [isj . 
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The main advantage of the Faddeev formulation is to allow a proper description of the different asymptotic states, 
each of them well separated by the corresponding Faddeev component and easily described in terms of their natural 
Jacobi coordinates. It is clear that the implementation of multichannel boundary conditions listed in (|8 UTT|) is 
impossible by means of one single function. Let us just however mention that this is not the only reason for using 
this approach, after all when working in momentum space there is no any need to implement any kind of boundary 
conditions. In fact there are also deeper reasons related to the mathematical structure of the underlying Lipmann- 
Schwinger like integral equations. The reader interested can take benefit in consulting [2^. [25j. 

Equations (I13|) allow to treat any physical non-relativistic quantum three-body problem with the only exception of 
break-up reaction with at least one attractive Coulomb interaction, despite numerous attempts [26| . 

Last but least, this formulation is also advantageous in the numerical calculations. Indeed, the Faddeev decom- 
position takes benefit of the symmetry properties of the system. As a consequence, t he Faddeev components have 
simpler structure than the total wave function itself and are therefore easier to interpolate numerically. 

The above described formalism has been extended to account for particles with spin, interacting via three-body 
and/or non local forces. Few years after their formulation, the Faddeev equations were generalized to an arbitrary 
number (N) of particles by Yakubovsky [2]. However their numerical solution is still limited to N=4 case and for 
the energies below the first three-body break-up threshold. In this manuscript we will not present the formalism of 
Yakubovsky equations and their numerical implementation techniques, interested reader may refer to f27| . Still in 
sections lIVlTVl we present some results, which were obtained by solving four-particle Yakubovsky equations. 

One should also mention that, although we have described Faddeev approach in coordinate space in terms of 
coupled differential equations, the equivalent formulation is available in momentum space using coupled integral 
equations. Both formulations have its own advantages, namely higher energy scattering is easier to treat in momentum 
space [l^ , while configuration space is more convenient for Coulomb related problems [3, [13, E^, HH] . 



III. CHARGED PARTICLE INTERACTION WITH H 



An interesting atomic physics problem is related to the polarization of neutral atoms. In the vicinity of a charged 
particle , the electronic density of a neutral atom A polarizes, giving rise to a central long range — A interaction 
which behaves asymptotically like 

Vpir) = -~ (14) 
2 

where ad is the dipole polarizability, characteristic of the atom (ad — 9/2 for H in the nip — > oo limit). This interaction 
is very shallow but if the charged particle is heavy enough to reduce the repulsive kinetic energy, it generates a rich 
spectrum of bound and resonant states. This can be qualitatively illustrated in a simple two-body approach by 
assuming an X^ — A two-body interaction of the form (ffH) with the replacement j29| : 

2 _2 / . 9 4 „ 3 27 2 27 27\ 
ad-^a{r) = ad - -e ( -|- -r"* + 9r3 + + + — j 

This particular choice of a(r), valid for scattering on H atom, regularizes (linearly) the polarization potential (jl4p at 
short distances, keeping also the right asymptotic behavior. 

We have plotted in Figure[5]the bound state spectrum of the tt — H system thus obtained. It contains a large number 
of states up to i = 6. In the case of a positively charged particle, X~^, these states would lay below the rearrangement 
threshold p + (e^X^). Figure [3] represents the fj, — H elastic cross section obtained in the same approximation, 
displaying very narrow resonances in high angular momentum states (l9| . 

In view of having reliable quantitative predictions, an exact calculation must be however performed, taking into 
account the internal structure of the atom. In the present Few-Body state of the art, this is only possible for the H 
and He atoms. The simplest case for which the calculations can be performed exactly is the scattering on atomic 
hydrogen. This constitutes a genuine three-body problem and a challenge for the Few-Body community. 

We have first considered the case of positively charged particle X^ with mx+ < mp, for they have at zero energy 
a single open channel corresponding to elastic scattering (see Figure S]) . The first exact three-body solutions of the 
X+(e~p"'") have been obtained in [l^]- We display in Figure [5] the X+ — H scattering length as a function of the 
mass of the projectile mx- The values corresponding to physical particles (/i,7r) are indicated by arrows. As one can 
see from this figure, the scattering length is divergent for some values of the incoming particle mass, indicating the 
appearance of an additional bound state in the {X,p, e) system. 
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FIG. 2: Bound state spectrum (atomic units) for tt-H 
in a two-body approach with a regularized polariza- 
tion potential (|14p . 



FIG. 3: Resonant fi~ H states in the same two-body 
approximation as in figure (2] 
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FIG. 4: DifTerent three-body reactions with atomic hydrogen and the existing particles. 



Of special interest is the scattering of protons. Not only for they constitute experimentally the most accessible 
physical particle but also because it corresponds, by chance, to one of the resonant cases displayed in Figure [S] for 
smaller projectile masses. 

Indeed we have computed the p — H scattering length for the S-wave symmetric (pp spin singlet) and antisymmetric 
(pp spin triplet) states. For the symmetric case we found = —29.3 a.u. while the antisymmetric one provided the 
value of at = 750 ± 5 a.u. The analysis of the nodal structure of the corresponding Faddeev amplitude indicated that 
such a big value is due to the existence of a first excited bound state with extremely small binding energy. By using 
the effective range expansion, its binding energy was found to be B=(1.135 ± 0.035) x 10^^ a.u. that is. ~ 30 neV. 

This state can be also viewed as the first excited vibrational level w = 1 of 2p(T„ symmetry in the molecular 
ion. A direct computation of this state using ad-hoc variational techniques [30j confirmed its existence and provided 



a more accurate value of the binding energy B — 1.085045 x 10 [31|. Further work showed that it is stable with 
respect to the relativistic and leading order QED corrections. Taking them into account, its binding energy is only 
slightly modified and becomes B = 1.082247 x 10"^ H. 

It is worth mentioning that this H2 first excited antisymmetric state exists also in the so called Landau p — H 
potential [s^. The latter consists in adding to the polarization term (I14p a repulsive one due to the Pauli principle 
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between protons in the spin triplet state and reads (in atomic units) 

^^(-)='?J?T-0 (15) 

The total potential Vl is regularized to a constant below Xc = 2.5. In its original formulation (ry = 1) it entails an 
excited S-wave state, although with binding energy two order of magnitude smaller than the exact three-body value 
and a p — H scattering length consequently larger. 
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FIG. 5: X'^ — H scattering length (in atomic units) as 
a function of the projectile mass (in MeV). The values 
corresponding to physical particles are indicated by 
arrows. 



FIG. 6: Cross section for the L=0 pH scattering in 
pp spin triplet state as a function of the energy (in 
atomic units). The body results (filled black square) 
are compared to those (red solid line) given by the 
two-body Landau potential [sj modified in order to 
reproduce the binding energy the first excited state. 



To our knowledge the Ht first vibrational 2ptT„ state above described constitutes the most weakly bound natural 
molecule ever predicted [59|. A direct computation provides a root mean squared radius R = 270 a.u. and its 
wavefunction has still sizeable values well beyond 1000 a.u.. That makes the state extremely unstable against any 
kind of perturbation. The low energy p — H scattering will be however totally dominated by the existence of this 
nearthreshold state. Corresponding S-wave elastic p — H cross section is shown in Figure [6l The Faddeev results 
(black filled squares) are compared to those (red solid line) obtained with the Landau potential ([T5|) . whose parameter 
was adjusted to 77 = 0.9254 in order to reproduce the exact 3-body binding energy. The huge values of the cross 
section are a consequence of the large triplet scattering length at . Notice that the resonant region is however limited 
to E < 10^^ a.u.. Below this energy range, a proton approaching an hydrogen atom will feel an object of nanoscopic 
size. 

Electron and positron scattering on Hydrogen atoms in their ground state presents one of the rare three charge 
structures for which experimental data are available. In fig. [7] we compare our calculated total scattering cross 
sections with the measured ones [s^, for energies below II(n = 2) threshold {E2 ~ 10.2 eV). One may see an 
excellent agreement between measured and calculated values. The only exception is in the e"*"-!! scattering at very low 
energies (E« leV), where the experimental result has a big error bar, probably due to the strong positron annihilation 
probability. 

We have displayed in figure [8] a zoom of the resonance region, close to the H(n = 2) threshold. One may identify 
three well-known narrow resonances in e~-H system, which are also confirmed experimentally [36l - l38| . One S-wave 
resonance is visible in e''"-II system. 

If the incoming particle is heavier than the proton, e.g — d, the rearrangement process 

X + (e-p) ^ p+{e-X+) 

is kinematically open and even at zero kinetic energy we are left with a two channel problem. By increasing the 
scattering energy, the number of open l-|-2 channels before the 3-body breakup threshold becomes infinite. 
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FIG. 7: e -H scattering cross section below H(n — 2) 
threshold: calculated values are compared with experi- 
mental ones jJSj. 



FIG. 8: A zoom in the resonance region of e^-H scatter- 
ing close to II(n = 2) threshold. 



The situation is even more interesting in the case of an incoming negatively charged which presents even at 
zero incoming energy a large numbers of open channel corresponding to Coulomb {X^ ,p) ground and excited states 
(see Figured]). All of them provides a very rich variety of phenomena, most of them remaining unexplored both from 
the theoretical as well as from the experimental point of view. 



IV. COLD MOLECULES 



An important simplification of the Faddeev equations may take place if one considers a system of identical particles. 
In this case, the three Faddeev equations, as well as the corresponding components, become identical and therefore 
the full solution of the problem is obtained by solving only one of them. Such symmetrized equation may be applied to 
study a system of identical bosons. This is the simplest few-body problem, nevertheless revealing genuine peculiarities 
of quantum many-body systems. 

Two of the striking phenomena manifested in few body physics are the Thomas collapse and the so called 
Efimov states [10, EH. 

The first of these two effects, the Thomas collapse, attests that an > 2 boson system will shrink into a state with 
an infinite binding energy, if the pairwise particle interaction is attractive and of zero-range. 

The second phenomenon, i.e. Efimov states, manifests in a series of weakly bound three-body states, which move 
below the two-body threshold by slightly increasing the two-body attractive interaction. These states can appear in a 
realistic physical system with finite range interaction, if one of its two-body subsystems is S-wave resonant, i.e. when 
the two body scattering length (ao) turns to be much larger than the interaction range {tq)- The interest in Efimov 
states has been fueled lately by the important progress of experimental low temperature physics and particularly by 
the recent discovery of an Efimov state in caesium atoms [42j and the subsequent, even more conclusive, experimental 
works of refs. [H-gl. 

Efimov physics involves systems having largely separated length (or momenta) scales and therefore turns to be a 
natural laboratory to test Effective field theories (EFT) . To analyze Efimov states it is then natural to introduce EFT 
expansions in terms of two small dimensionless parameters ro/ao and fcrg, k being a characteristic center-of-mass 
(cm.) momenta. 

It is a common believe that, at the leading order of EFT, all the low-energy properties of a three-boson system are 
set by two parameters: one two-body parameter - like two-boson (dimer) binding energy or scattering len gth - and 
one three-body parameter - like three-boson (trimer) binding energy or particle-dimer scattering length [47l l48j. This 
means that whenever in a three-body system the particle intcrdistance R satisfies the Efimov condition, vq <^ R <^ \ao\, 
its low-energy properties must be interaction independent and be governed by the only two low energy aforementioned 
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FIG. 9: Evolution of the three- and four-boson binding energies with the Gaussian cutoff A. 

parameters. 

More recently, it has even been demonstrated that the properties of the four-boson system (tetramer) are determined 
by the same two parameters and no additional four-body scale is required to establish universal relations between 
three- and four-boson observables [49, 50J. These universality studies have been also extended to fermionic three- and 
four-body systems, recovering series of phenomenologically observed correlation rules in multi-particle system [sij . 
Most of these relations have been derived relying on purely attractive (and often contact) two-particle interactions, 
while N > 2 systems are balanced by introducing repulsive three-particle force to prevent Thomas collapse . 

We present in this section a slightly different quantum-mechanical approach, based on a family of short range 
separable, rank-1 and rank-2, interactions having the same two-boson scattering length. 

These potentials have been derived somehow mimicking interaction between "^He atoms, i.e., choosing the boson 
mass h? /m — 12.12 K-A^ and fixing the two-boson scattering length to ao =104.0 A. Nevertheless one can not pretend 
with these results a realistic description of *He multimers, since these separable potentials do not retain well known 
properties of the effective interaction of two inert neutral atoms, like strong repulsion at the origin and weak van-der 
Walls attraction in the asymptote. 

The family of two-boson separable (and thereby nonlocal) potentials we use have the form: 

V = ^\g^)kj{9J\ (16) 
y 

where n,. is the rank of the considered potential. For simplicity, these potentials are restricted to S'-wave only. In close 
analogy with most of the EFT potentials [i^, the form factors \gi) are chosen as Gaussian. In the configuration-space 
they are given by 

(r|5.) = (A,/y2)3e-(^-/2)^ (17) 

The rank-1 potential is defined by two parameters: An and Ai. To preserve the two-body scattering length at 
oo =104.0 A, these two parameters must satisfy An = [7rm(l/ao — Ai/-\/27r)/2]~^. 

The rank-2 potential is obtained by choosing Aj = jA with A in the range [0.133,0.305] A-i, i.e. 13.8 < gqA < 31.7, 
and determining the remaining three strength parameters Ay (we assume A12 = A21) by a fit of the calculated 
observables (with a typical four digit accuracy): dimer binding energy B2 =1.318 mK, scattering length oq =104.0 
A and two-particle scattering phase-shifts up to about 50 B2 cm. energy. The broad scattering energy interval 
included in the fit, roughly corresponding to the natural energy scale fi^/mrg [4^, guarantees that in the effective 
range expansion not only but also higher order terms are well retained. Note that the phaseshift set we have used 
to fix the parameters of the rank-2 potential differs from the one obtained by realistic interaction (for instance our 
effective range vq = 15.0 Ais about two times larger). 

The dependence of the binding energies of L=0+ three- and four-boson states as a function of the Gaussian cut-off 
A fEd] is displayed in fig. [S] for both rank-1 and rank-2 potentials. By binding energy i?„ of a n-particle multimer we 
denote the difference between the sum of the n free constituent particles masses and the total multimer mass. 
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The rank-1 results are consistent with the EFT predictions, based on contact 2- and 3-body interactions. There 
exist two trimers, the ground state which is strongly bound (its binding energy exceeds the dimer one by more than 
a factor 10), and the first excitation which remains close to the dimer threshold. There are also two tetramers, 
the ground state being strongly bound relative to trimer one and the excited tetramer which consist in one boson 
weakly attached to the trimer ground state. The binding energies of these trimer and tetramers reveal almost a linear 
correlation pattern for rank-1 potential, see Fig. [101 and are in agreement with the EFT results of ref. (5ll |. 

The situation is much more complicated when using the rank-2 potential. The binding energies B^^^ / B2 are slowly 
decreasing functions of the cut-off A in the range oqA < 26.0. In this region, all potentials support one ground and 
one excited state for the trimer as well as for the tetramer, in agreement with the EFT approach. Moreover the trimer 
and tetramer binding energies follow closely the correlation curves of rank-1 values. However, for larger A values, 
the rank-2 potential results show a nontrivial behavior despite the fact that all the potential parameters \ij depend 
smoothly on A. 

First, at oqA = 26.5, the binding energy of the excited tetramer i?4^'' detaches from b'"^\ It increases drastically 
and stabilizes for a moment when it almost attains the tetramer ground state binding energy B^^^ that, with a tiny 
delay, also starts to increase rapidly. At the same time, close to oqA = 26.75, the second excited tetramer state 
appears with a binding energy i?^^-* only slightly larger than Bg"-* . In this region of harsh tetramer variations, the 
trimer properties remain practically unchanged. 

The trimers exhibit qualitatively similar behavior, although at larger A values, around agA = 27.5. The binding 

energies of the trimer ground and first excited state, Sg"' and B^^\ start to grow very rapidly; the former increases 
much like Bf^ did while the latter stabilizes for a while. As b'"^^ increases, the second excited tetramer state disappears 
below the ground trimer threshold at oqA = 28.14, whereas the first excited tetramer stalls for a while letting the 
trimer threshold B'f^ to approach again. Furthermore, at oqA = 28.12 the second excited trimer state appears. 
Unlike the strongly bound trimers, this state behaves like an Efimov one, i.e., it slides under the dimer threshold if 
the two-boson interaction is made stronger. The second excited trimer loses that property above oqA = 30.7 where its 
binding energy b'^\ much like b'^^ at oqA = 27.5, starts a phase of a rapid increase. However, the Efimov property 
is recovered by the third excited trimer, appearing at oqA — 31.7. It is interesting to note that the second excited 
tetramer shows a similar behavior, i.e. it also disappears if the two-boson interaction is made stronger. 

While rank-2 potential results seems to be rather confusing, the explanation of such a nontrivial behavior of a few- 
boson system is rather simple. For this aim one should simply understand the A-evolution of the nonlocal interaction. 
The shape of = 1 potential is very simple, and can be represented by a single attractive bell-shape well in (r, r') 
plane, the width of the bell being constrained by the size of the cutoff A. For small A values, the rir = 2 potential has 
also single and rather flat attractive plateau in the (r, r') plane. However for larger A values, this plateau splits into a 
rather complicated saddle-like surface with two asymmetric attractive regions separated by two symmetric repulsive 
regions. A narrow attractive region is formed at the origin and it deepens when increasing A, whereas the second, 
distant, attractive region is much wider and flatter. Due to the large kinetic energies required to squeeze the particles 
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into the narrow region at the origin, it is energeticahy preferable for the multi-boson system to stay in the distant 
attractive plateau. At this moment, the rij. — 2 potential presents qualitatively the same features as = 1, having 
a single attractive region. However, for some values close to aqA = 26.0, the attractive region at the origin becomes 
deep enough to accommodate the four-boson system, and some kind of a phase transition occurs. The trimer has 
only three pair interactions compared to six for the tetramer. Therefore the corresponding phase transition in the 
three-boson system takes place at larger A values, when the attractive region near the origin is even more deepened. 

One can probably question the reality of such "exotic" nonlocal potentials with several attractive regions. Nev- 
ertheless, one must admit that some realistic nucleon-nucleon potentials [52, 53] have even more complex structure 
than the one exhibited in our model. Effective interaction between two complex molecules can however reveal even 
much more complicated pattern. 

The presence of finite range (rg > 0) interactions does not guarantee that the Thomas-like collapse will not occur 
in multi-boson system. The interaction can have a short-range or an off-shell structure that may be ignored with the 
low-energy probe. Such structure can have no or little effect for low-energy observables of an N particles system, 
but may be exploited in the N -\- 1 particle one which, having more interacting pairs, can compensate larger kinetic 
energies and thus regroup itself into a shorter range domain. This fact complicates the possibility to make universal 
predictions about multi-particle system, in particular about its deeper lying states, based only on the low energy 
properties of its subsystems, unless one is sure that the eventual size of the few-body system satisfies the Efimov 
condition i? >> rg. Of course, this condition is never satisfied a priori, although one may expect that it must be valid 
for higher order excited states fs^j. 



V. *HE MULTIMERS 



In this section we turn our attention to a realistic description of small structures of ^He atoms. The two electrons 
of the He atom close the Is shell, determining its spherical symmetry as well as its inert chemical properties. Two He 
atoms repel strongly when approaching to each other, but a very weak Van der Vaals attraction is exhibited at large 
distance, resulting in a shallow attractive pocket with a maximal depth of RilO.9 K, centered at Rue-He ~ 3 A. 

This weak attraction is responsible for the fact that at very low temperatures both bosonic '''He but also fermionic 
■^He become liquid. It is also responsible for the existence of a loosely bound "^He dimer, with binding energies of 
B fa 1 mK and sizeable "^He- "^He scattering length ao ~ 100 A, which has been experimentally observed few years 
ago [Hill. 

When trying to study theoretically N>2 "^He clusters, one is faced to a serious theoretical problem. The standard 
boundary conditions in the Dirichlet form (|8]), which imposes to the FY components to vanish at the origin, turns 
to be impractical when considering structures of He atoms. The strong hard-core repulsion describing the inner 
part of the He-He potential, which dominates the interaction between two He atoms at relative distance Rne-He ~2 
A, introduce severe numerical complications. The relevant attractive matrix elements which are responsible for the 
binding fade away in front of these huge repulsive hard-core terms, thus causing numerical instabilities. On the other 
hand, such strong repulsion close to the origin physically, simply reflects the fact that the two He atoms cannot get 
close to each other, i.e. closer than some relative distance r = c inside the core. 

The wave function of three He atom system, for instance, should vanish in the region of three-particle space M^, 
which is the interior of the three multidimensional surfaces Xi = c, where Xi is the distance between the particle 
j and k defined in ([T]). Therefore one way to attain numerical stability would be to neglect the solution in the 
strong repulsion region. In practice, as it has been shown by Motovilov and Merkuriev (STj . the presence of infinitely 
repulsive interaction at Xi < c, can be formulated in terms of boundary conditions for the Faddeev components, which 
are implemented by setting: 

[E -Ho~ V{xi)] ^^{x^,y^) = for X, < c, 

^^{x^,y^)=0 for .T, = C. (18) 

A similar approach may be applied also to condition Faddeev- Yakubovsky equations for four-body system, see (STj . 

We compare in Table |T] the calculated properties of three and four "^He atom structures. The binding energies 
of the ground and excited states (in mK) together with atom-dimer and atom-trimer scattering lengths (in A) are 
listed. The realistic model is based on the interaction between two He atoms developed by Aziz and Slaman [ssf . 
popularly referred to as LM2M2 potential. These realistic calculations are compared to the values obtained using the 
rank-1 potential from the previous section, whose parameters (A, A) have been adjusted in order to reproduce the 
dimer and ground state trimer binding energies obtained with the LM2M2 potential. One may see that the rank-1 
approximation reproduces well the excited trimer binding energy Bg^'' and the atom-dimer scattering length Oq^ ^\ 
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TABLE I: Predictions for the ''He multimer properties. Binding energies (in mK) and scattering lengths (in A) obtained using 
realistic LM2M2 interaction and rank-1 separable model of section ITVl 



Pot. 


5(0) 


^(1) 

-°3 


5(0) 


o(l)_o(0) 
^4 ^3 


{2-1} 
% 


{3-1} 


LM2M2 


126.39 


2.2680 


557.7 


1.087 


115.56 


103.7 


rank-1 


126.39 


2.2758 


597.9 


3.16 


114.77 


67.38 



However the description of the four atom system deteriorates considerably. Probably the main reason is the difference 
of the off-energy shell properties between these two interaction models: a purely attractive rank-1 potential and the 
LM2M2 one containing a strong hard core region. The four boson system, being a more compact structure, test more 
strongly the hard core region than the three body one, which is also reflected by the fact that the rank-1 potential 
provides more binding than the LM2M2 one. 

It is worth mentioning that a direct calculation of the ^He tetramer excited state still represents a challenging 
numerical task. This state is very weakly bound and its wave function extend over several hundreds of A. One is 
therefore forced to use a very large and dense grid, ensuring at the same time enough accuracy to trace the small 
binding energy difference with respect to the trimer ground state. Nevertheless the vicinity of the tetramers excited 
state to the ^He- ^Hea continuum makes possible the extraction of its binding energy from the low-energy scattering 
results as explained in [l^ . 

VI. CONCLUSION 

The Faddeev-Yakubovsky equations provide a rigorous quantum mechanical formulation of the non relativistic few 
particle problem, taking into account all the possible asymptotic states of the system. First formulated in momentum 
space and for short range pair-wise interactions, they have been generalized to accommodate three-body forces as well 
as long range Coulomb potentials. 

Aside from its well defined mathematical structure, the Faddeev-Yakubovsky formalism turned to be very useful 
in numerical calculations. A large variety of problems in hadronic, nuclear, atomic and molecular physics have been 
solved since its appearance. Still limited to N=3 and 4 systems, they can be in principle extended to an arbitrary 
number of particles. They provide very accurate results for the bound state problems but the genuine feature of this 
formalism lies on its ability to treat scattering and bound states on the same footing, thus enabling a complete and 
consistent analysis of the few-body system under consideration. 

The configuration space formulation of Faddeev-Yakubovsky equations, made possible once the boundary conditions 
are well established, is the best adapted to deal with atomic physics problems. Some selected examples of them have 
been presented in this contribution covering aspects of cold atomic and molecular physics. We have considered 
the scattering of positively charged particle on atomic hydrogen, systems of cold atomic molecules and bound and 
scattering states of *He atoms. In particular, we would like to emphasize the result concerning the p — H scattering: 
the huge scattering length in the proton-proton spin triplet state at ~ 750 atomic units. 

Even by restricting ourselves to the 3- and 4-body cases, a large variety of unsolved problems remains to be 
considered. Among them we would like to mention: 

• The Coulomb reactions including break-up and/or many open channels, like those presented in Figure |4] 

• The many interesting antiproton physics processes, like p + d, pH — > + (pp)* , p + (e^e^) — > _ff + e~ 

• Four-body Coulomb problems, like first principle calculations of positron-positron scattering 

• For-Body break-up reactions involving 1+2+1 and 1+1+1+1 particle channels 

They constitute altogether a very rich and challenging program for the coming years, both from the theoretical and 
experimental physics point of view. 
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